Dynamics of semi-flexible polymer solutions in the highly entangled regime 
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We present experimental evidence that the effective medium approximation (EM A), developed by 
D.C. Morse [Phys. Rev. E 63, 031502, (2001)], provides the correct scaling law of the macroscopic 
plateau modulus G" oc p'^^^Lp^^^ (where p is the contour length per unit volume and Lp is the 
persistence length) of semi-flexible polymer solutions, in the highly entangled concentration regime. 
Competing theories, including a self-consistent binary collision approximation (BCA), have instead 
predicted G" oc p''^^ Lp^^^ . We have tested both the EMA and BCA scaling predictions using 
actin filament (F-actin) solutions which permit experimental control of Lp independently of other 
parameters. A combination of passive video particle tracking microrheology and dynamic light 
scattering yields independent measurements of the elastic modulus G and Lp respectively. Thus we 
can distinguish between the two proposed laws, in contrast to previous experimental studies, which 
focus on the (less discriminating) concentration functionality of G. 
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Despite their importance to soft-matter physics, biol- 
ogy and industrial processing, the viscoelastic properties 
of semi-flexible polymer solutions are still not well un- 
derstood and a basic analytical model has not yet been 
agreed upon. All current models describing the viscoelas- 
tic properties of semi-flexible polymer solutions are elab- 
orations on the early models of Doi and Edwards [3, . 
They developed two full theories of the entangled state 
for two extreme cases: completely flexible 0] and rigid- 
rod 3 polymers. Solutions of semi-flexible polymer, 
that lie between those extremes, have many regimes of 
viscoelastic behavior (requi ring many theoretical mod- 
els in, i, S i, 0, S i, [3, [ii 0), depending on the 
polymers' degree of rigidity (described in terms of persis- 
tence length Lp), on their contour length L, and on the 
concentration (from dilute to highly entangled regimes). 
We shall focus on highly entangled isotropic solutions of 
semi- flexible polymers, with L/ Lp ~ 1. In particular, we 
study the range of concentration (« 0.1-1 mg/mL in this 
case) where the geometrical mesh size Lm is much less 
than Lp, and the tube diameter and entanglement length 
are also expected to be much less than Lp. This range 
of concentration was deflned by Morse [l| as the tightly 
entangled regime, and is particularly relevant to many 
biological and industrial polymeric fluids. In order to de- 
scribe the viscoelastic behavior of the polymer network in 
this range of concentrations, Morse developed two ana- 
lytical approximations describing the conflnement forces 
acting on a randomly chosen test chain embedded in a 
"thicket" of uncrossable chains: the binary collision ap- 
proximation (BCA) and effective medium approximation 
(EMA). In fact, the scaling relation resulting from the 
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BCA had previously been obtained by several others au- 
thors but Morse has also estimated the pref- 
actors. So, prior to the introduction of the EMA, there 
was broad agreement regarding the scaling law. The ap- 
proximations are summarized as follows. 

(i) The binary collision approximation gives a 
rather detailed description of the interaction of a test 
chain with individual nearby medium chains, but neglects 
any effects arising from the collective elastic relaxation 
of the network. It yields the following expression for the 
elastic modulus: 

G w O.AQkBTp'^'^Lp'^'^. (1) 

(ii) The eflfective medium approximation starts 
from a very different point of view, by treating the net- 
work surrounding the test chain as an elastic continuum 
with a shear modulus equal to the self-consistently deter- 
mined plateau modulus of the solution, and the test chain 
as a thread embedded in this medium. The expression 
thus obtained is 

G « 0.82A:Br//3x-i/3 (2) 

Comparison of the above scaling predictions raises the 
question of which theoretical approach (if either) bet- 
ter describes the viscoelastic behaviour of semi-flexible 
polymer solutions in the tightly entangled concentration 
regime. Existing experimental measurements of the con- 
centration dependence of the plateau modulus are 
not sufflcient to answer that question, because the two 
putative values of the exponents are numerically quite 
close (G oc p^-^ vs. G oc p^'^^), so that they both fit the 
experimental data with reasonable accuracy [ij . 

To test the scaling predictions, we experimentally anal- 
ysed solutions of actin filaments (F-actin), a semi-fiexible 
polymer derived from muscle tissue. As we shall show, F- 
actin has the useful property that its persistence length 
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can be controlled by varying only the ionic properties 
of its solvent, without altering other system parameters 
such as solvent viscosity, polymer concentration or molec- 
ular weight, thus allowing us to discriminate between the 
EMA and BCA models. Hence, this is a rare example of 
biology helping to answer questions of interest to physics, 
rather than vice versa. 

At low ionic strength in vitro, actin exists in the 
monomeric (globular) G-actin form. G-actin is roughly 
spherical with a diameter of about 5 nm. When the ionic 
strength of a G-actin solution is increased to a physiologi- 
cal value (0.1 M), G-actin self-associates, to form F-actin, 
which is characterized by a persistence length of 2-20 ^m 
(depending on the buffer used) [13, and a diameter of 
approximately 8 nm [10, Tl- In order to produce actin 
filaments with different mechanical properties (Lp), two 
sets of F-actin solutions (named System 1 and 2) were 
prepared using two different buffer recipes, which will be 
referred to hereafter as F-bujfer 1 and F-huffer 2. In both 
the cases, the initial solutions of G-actin were prepared 
with the same buffer recipe {G-buffer). Polymerization 
was initialized by increasing the ionic strength of the G- 
buffer to that of the F-bujfers. The buffer recipes are as 
follows: G-buffer: 0.2 mM ATP, 0.2 mM CaC12, 2 mM 
Tris-HCl, 0.5 mM DTT, pH 8.0. F-buffer 1 : 50 mM KCl, 
2.0 mM free MgC12, 5 mM Tris-HCl, 1 mM ATP buffer 
pH 7.5. F-buffer 2: 25 mM KCl, 1.0 mM free MgC12, 1.0 
mM EGTA, 10 mM MOPS buffer pH 7.0. 

In order to obtain independent measurements of the 
complex modulus G* and of Lp, the two sets of F-actin 
solutions were each investigated by two different tech- 
niques: passive video particle tracking microrheology 
(PVPTM) and dynamic hght scattering (DLS). 

The first technique, PVPTM, exploits the relationship 
between the viscoelastic properties of a fluid under in- 
vestigation, and the mean-square displacement (MSD) of 
probe particles (of radius a), suspended in the fluid and 
executing Brownian motion (Fig.[T]). This relationship is 
given by the Generalized Stokes-Einstein relation. 



G{s) = kBT/T:as{Af\s)), 



(3) 



where (Af^(s)) and G{s) are the Laplace-transforms of, 
respectively, the mean square displacement and the time 
derivative of the shear modulus. Moving from Laplace 
space to Fourier space is accomplished by substitut- 
ing the Laplace frequency s with iuj, so that the real 
and imaginary parts of G{iuj) correspond to the storage 
and loss moduli, G'{lo) and G"{uo) respectively, with u) 
the Fourier frequency (Fig. [2]). As probe particles, we 
used carboxylate-coated polystyrene beads of diameter 
0.489 /xm. It is well known that microrheology can fail 
to emulate macroscopic results if the probe particle's size 
and surface chemistry are incorrectly chosen. Microrhe- 
ological measurements of viscoelastic moduli are often 
found to match the scaling laws found by macroscopic 
rheology, whilst disagreeing by a constant factor in ab- 
solute magnitude. In the present case, we only require 
reliable measurements of scaling exponents, not of abso- 
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FIG. 1: Mean-square displacement (MSD) versus lag-time, 
for 0.489 /im-diameter beads in F-actin solutions at different 
concentrations [mg/mL] (System 2). The solid line represents 
the MSD of the beads suspended in water at 25°C. 
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FIG. 2: Storage and loss moduli versus frequency, for F-actin 
solution at a concentration 1 mg/mL (System 2). The lines 
are guides for the eye. The same features are also exhibited 
by System 1 



lute values. Nevertheless, Fig. [3] demonstrates that the 
absolute values of our measurements are at least as ac- 
curate as macroscopic rheometric data in the literature, 
so that we can be confident of their validity. 

In all our PVPTM measurements, we found very good 
agreement with Morse's scaling predictions, for both 
G'{uj,c) and G"(liJ,c) as functions of frequency and con- 
centration. Figure [3] shows the plateau modulus G^ (de- 
termined as the value of the storage modulus G'(w) at the 
frequency for which the ratio G" {lo) / G' {uj) is minimum) 
versus F-actin concentration, for Systems 1 and 2 in the 
tightly entangled regime (w 0.1-2 mg/mL ). It is clear 
from the figure that, as previously found [l|, the con- 
centration dependence of G° is not able to discriminate 
between the two predictions (EMA and BCA). However, 
using the microrheology data, the two theories yield very 
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FIG. 3: Plateau modulus G" versus F-actin concentration, 
for Systems 1 (circle) and 2 (square). The open triangles 
are bulk rheology results taken from ref. [l^. The short 
dash and dash-dot lines are linear fits to the data on the 
log- log plots (power laws); the dot and dash lines show the 
BCA (G° oc c'^^) and EMA {G° cx c^^^) scaling predictions 
respectively. The inset (System 2) shows the actual existence 
of a different power-law regime at concentrations lower than 
the tightly entangled regime (c.f. Ref. [T^l. 

different predictions for the ratio of persistence lengths 
in the two systems: 



BCA Lpi/Lp2 = {G°2/Glf = 18 ± 3 



(4) 



EMA Lpi/Lp2 = {G°2/Glf ^ 6.2 ± 0.7 (5) 

Hence, we can distinguish between the two models by 
measuring the actual persistence lengths. For that pur- 
pose, DLS measurements were performed on some of the 
samples. 

Kroy and Prey [l7| provide a simple analytical ex- 
pression for the time-dependent decay of the dynamic 
structure factor g*^^^ {q, t) of semi-flexible polymers in 
semi-dilute solution, where q — [(47m)/Ao] sin(6'/2) is 
the magnitude of the scattering wave-vector, defined as 
the difference between the incident and scattered wave- 
vectors, n is the refractive index of the solvent, Ao is the 
wavelength of the laser in vacuo, and is the scatter- 
ing angle. They showed that, at sufficiently long times 
{t » (gLp)~^/^7~^), the dynamic structure factor re- 
duces to a simple stretched exponential, 

g(i)(g,t) =g(i)(g,0)exp [~T{l/i){^,tf^y37r] , (6) 

where r(l/4) = 3.62561 and decay rate jq given by 

7, = fcsTq^/^ [5/6 - Hqa^)] /^TrVsL^^^ (7) 

where 77^ is the solvent viscosity. They also calculated the 
form of the dynamic structure factor in the limit f — > 0, 
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FIG. 4: Double logarithm of the normalized dynamic struc- 
ture factor g^^\q,t)/g^'^\q,0) versus logarithm of time, at a 
scattering angle of 90°, for a solution at F-actin concentra- 
tion c = 0.4 mg/mL (System-1). Similar data are obtained 
for System 2. The lines are guides for the eye. 



which can be used to measure the microscopic hydrody- 
namic lateral diameter ah of the polymer chain: 



/'\q,t)=g^'\q,0)e^p{-j't) 



with initial decay rate 



yO = fcsTq^ [5/6 



ln(ga,J] /Qn'^ris 



(8) 



(9) 



The validity of Eqs. ([6]-[8]) requires [l3| the conditions 
ah << q~^ << Lp ^ L and q~^ << Lm, all of which 
are satisfied by all of our samples, since ~ 9 nm, 
« 125 — 33 nm (where only the high-g data will be 
used in the final result), while Lp and L are of order of 
microns and L„i ~ 0.51 — 0.16/im [T]| . It is clear from 
Fig. d] that Kroy and Frey's calculations [l3| are appli- 
cable to our experimental system, despite the scatter of 
data at early times. Indeed, 5^^-' (q, t) decays initially as a 
simple exponential (identified by the line with unit slope 
in Fig. 2]) and exhibits stretched exponential behavior at 
longer times (line with slope 3/4). 

The hydrodynamic diameter ah was obtained from 
the data by evaluating the initial decay rate, 7'^, via a 
linear fit to the logarithm of the normalized dynamic 
structure factor, in a time window between 10~^ s and 
4 X 10"^ s. From Eq. ^ the non-dimensionalized ini- 
tial decay rate, r° — Q-K^rfs^i^ / {ksTq^), was then fitted 
by [5/6 — ln(ga/i)] (Fig. [5]) to obtain ah- For all the F- 
actin solutions investigated (as with other large polymers 
[ill), we found very good agreement between the indirect 
measurement of the lateral hydrodynamic diameter (av- 
eraged value Oft" = 9 ± 3 nm) and that measured directly 
by transmission electron microscopy (TEM) image anal- 
ysis [oh = 9.2 ±0.1 nm; see Fig. [S] insert) . 
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FIG. 5: Normalized initial decay rate T^' versus scatter- 
ing wave-vector q, and the best fit of Eq. Q using ah as 
a free parameter, for a solution at F-actin concentration 
c — 0.1 mg/mL (System-2). Insert: Diameter distribution 
from TEM image analysis. 
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FIG. 6: Persistence length Lp versus wave-vector q for a solu- 
tion at F-actin concentration c — 0.2 mg/mL (System-1) (de- 
rived using ah = 9.2 nm). 



With the hydrodynamic diameter now determined, it 
can be used in Eqs. ^ and (O to estimate the persis- 
tence length Lp. We adopt an average value = 9.2 nm 
for both systems. Equations ([6]) and ([7]) were fitted to the 



measured time-dependence of the dynamic structure fac- 
tor, for various scattering vectors q, with Lp as the only 
fitting parameter. The results for System 1 are shown in 
Fig. [HI One might expect the result to be independent 
of q but, as noted in Ref. the equations apply only 
to scattering vectors large enough to resolve single fila- 
ments. For scattering vectors smaller than the inverse 
mesh size Lm, the structure factor is averaged over sev- 
eral filaments and Eqs. ([6]) and (O no longer hold. It is 
therefore the large-g asymptote that marks the true per- 
sistence length in Fig. [B] We thus consider only averaged 
values of Lp for q >22 /im~^, as reported for the two sys- 
tems in Table U where three independent measurement 
of the ratio of persistence lengths are all in agreement. 

In conclusion, the true ratio of persistence lengths in 
the two systems is consistent with the value predicted 
by Morse's EMA approximation (equation [2]) 1!}, but is 
more than four standard deviations away from that pre- 
dicted by the more established BCA scaling law (equa- 
tion [T|) . We have thus resolved the controversy, and es- 
tablished that the effective medium approximation more 
accurately models the tightly-entangled regime of semi- 
flexible polymer solutions. The strength of this result 
rests on the clear agreement between two completely in- 
dependent experimental techniques. 
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Actin concentration 


Lpi 


Lp2 


Lpi /Lp2 


(mg/mL) 


(/im) 






0.1 


10 ± 1 


1.88 ±0.05 


5.3 ±0.6 


0.2 


12.2 ± 0.8 


2.24 ±0.15 


5.4 ±0.5 


0.4 


18 ±2 


3.16 ±0.09 


5.7 ±0.7 



TABLE I: Persistence lengths in Systems 1 and 2. We also 
found a weak concentration dependence of Lp oc c^/'^*""* 
for both Systems. This functionality is consistent with the 
definition of Lp, which is expected to grow from Lp — hi/ksT, 
in dilute solution (where k is the bending modulus), up to 
Lp = oo, in the nematic phase. 
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